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ABSTRACT 

Having found that type 1 Seyfert nuclei have excess [Fe VII] A6087 emission 
with respect to type 2s, Murayama & Taniguchi have proposed that the 
high-ionization nuclear emission-line region (HINER) traced by the [Fe VII] 
A6087 emission resides in the inner wall of dusty tori. The covering factor of 
the torus is usually large (e.g., ~ 0.9). Further, electron density in the tori 
(e.g., ~ 10^^® cm~^) is considered to be higher significantly than that (e.g., 
~ 10^^^ cm~^) in the narrow-line region (NLR). Therefore it is expected that 
the torus emission contributes to the majority of the higher-ionization emission 
lines. Taking this HINER component into account, we have constructed new 
dual-component (i.e., a typical NLR with a HINER torus) photoionization 
models. Comparison of our model with the observations show that that if the 
torus emission contributes ~ 10 % of the NLR emission, our dual- component 
model can explain the observed high [Fe VII] A6087/[O III] A5007 intensity 
ratios of the Seyfert Is without invoking any unusual assumptions (e.g., the 
overabundance of iron). 

Subject headings: galaxies: active - galaxies: Seyfert - galaxies: emission lines 

1. INTRODUCTION 

It has been often considered that emission-line regions around active galactic nuclei 
(AGN) are photoionized by the non-thermal continuum radiation from central engines 
(Davidson 1977; Kwan & Krolik 1981; see for a review Davidson & Netzer 1979; Osterbrock 
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1989). However, this photoionization scenario has been sometimes confronted with 
several serious problems; 1) any single-cloud photoionization models underpredict the 
[O III] A4363/[0 III] A5007 intensity ratio (e.g., Filippenko & Halpern 1984), 2) any 
(ionization-bounded) single-cloud models cannot explain the large temperature difference 
between T[oiii] and Tjnh] (e.g., Wilson, Binette, & Storchi-Bergmann 1997), 3) any 
single-cloud photoionization models underpredict higher ionization emission lines including 
ultraviolet emission hues such as C III] A1909 and C IV A1549 (e.g., Pelat, AUoin, & 
Bica 1987; see also Dopita et al. 1997), and 4) any single-cloud photoionization models 
require the overabundance in some particular elements such as nitrogen and iron (e.g.. 
Ho, Shields, & Fihppenko 1993; Hamann & Ferland 1992, 1993). In order to reconcile 
the above dilemmas, multi-component cloud models have been proposed (Stasihska 1984a, 
1984b; Filippenko & Halpern 1984; Binette 1985; Ferland & Ostcrbrock 1986; Binette, 
Robinson, & Couvoisier 1988; Binette et al. 1996; Komossa & Schulz 1997; Ferguson, 
Korista, & Ferland 1997a; Ferguson ct al. 1997b). Although this approach seems to be the 
natural extension of the single-cloud photoionization models, there is no guiding principle 
for the multi-component modeling because we do not have any unambiguous observational 
constraint on the geometrical structure of the ionized clouds in AGN. 

Recently Murayama & Taniguchi (1998, hereafter MT98) have found that type 1 
Seyfert nuclei (Sis) have excess [Fe VII] A6087 emission with respect to type 2s (S2s). 
The Sis exhibit broad emission hues which are attributed to ionized gas very close to the 
central engine (e.g., ~ 0.01 pc), whereas the S2s do not show such broad lines. The current 
unified model of AGN explains this difference as that the broad line region in the S2s is 
hidden from the line of sight by a dusty torus if we observe it from a nearly edge-on view 
toward the torus. Therefore, the finding of MT98 implies that the high-ionization nuclear 
emission-line region (HINER: Binette 1985; Murayama, Taniguchi, & Iwasawa 1998) traced 
by the [Fc VII] A6087 emission resides in the inner wall of such dusty tori. Since the covering 
factor of the torus is usually large (e.g., ~ 0.9), and the electron density in the tori (e.g., 
~ 10^~^ cm~^) is considered to be significantly higher than that (e.g., ^ 10^~^ cm~^) of the 
narrow- line region (NLR), the contribution from the torus dominates the emission of the 
higher- ionization lines (Pier & Voit 1995). Taking this HINER component into account, we 
construct new dual-component (i.e., a typical NLR with a HINER torus) photoionization 
models and examine whether these models can explain the observations without invoking 
any unusual assumptions. 



2. PHOTOIONIZATION MODELS 
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2.1. Single Cloud Models 

In order to investigate the HINER of Seyfert galaxies within the framework of 
photoionization, we first perform single-cloud photoionization model calculations using the 
spectral synthesis code CLOUDY version 90.04 (Ferland 1996) which solves the equations 
of statistical and thermal equilibrium and produces a self-consistent model of the run of 
temperature as a function of depth into the nebula. Here we assume that a uniform density, 
dust-free gas cloud with plane-parallel geometry is ionized by a power-law continuum 
source. 

The parameters for the calculations are 1) the hydrogen density of the cloud (nn), 
2) the ionization parameter ([/), 3) the spectral index of ionizing continuum between 10 
//m and 250 A (a; /j, oc and 4) the chemical compositions. We adopt the following 
continuum spectrum taking account of the actually observed one in AGNs (cf. Ho, Shields, 
& Filippenko 1993); 1) a = 2.5 for A > lO/xm, 2) a = -2 - -1 for 250 A< A < 10 /xm, 
and 3) a = — 2 for A < 250 A. The total amount of metals is expressed as Z in units of 
Zq. In this Letter, the abundance of each element is defined as a number fraction relative 
to hydrogen and is expressed by bracketed quantities ; e.g., [Fe] = n(Fe)/n(H). All the 
abundances are in units of solar values. However, taking into account the many lines of 
evidence in favor of a nitrogen overabundance (cf. Storchi-Bergman et al. 1996; Hamann 
& Ferland 1992, 1993), we adopt twice the solar nitrogen abundance as the standard 
value in our model calculations following Ho, Shields, & Fihppenko (1993). Namely, all 
elements have solar values except for nitrogen whose abundance will be twice the solar 
value. Adopted solar abundances relative to hydrogen are taken from Grevesse & Anders 
(1989). The calculations were stopped when the temperature fell to 4000 K, below which it 
was verified that little optical emission took place. 

We are not going to present detailed description of these single-cloud models because 
the main purpose of this Letter is to present new dual-component photoionization models. 
Therefore, we give a summary of the single-cloud models briefiy (see for details Murayama 
1998) and describe the difficulty of the models. We have performed several model runs 
covering the following parameter ranges: ei) a — —2, —1.5, and —1, b) logC/ = —2.5, —2, 
-1.5, and -1, c) lognn (cm'^) = 4, 5, 6, and 7, d) [Z] = 1, 2, 5, and 10, and e) [Fe] = 1, 
2, 5, and 10 keeping the other elemental abundances to the solar values. These parameter 
ranges were adopted in order to reproduce the observational properties of the HINER. We 
have confirmed that the intensity of [Ca V] A6087 is always less than 5% of that of [Fe VII] 
A6087 in all the cases, providing the validity of the use of [Fe VII] A6087 emission line. 

Using several newly defined diagnostic diagrams (e.g., [Fe VII] A6087/[O III] A5007 vs. 
[O I] A6300/[O III] A5007, [Fe X] A6374/[Ne V] A3426 vs. [Fe VII] A6087/[Ne V] A3426, and 
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[Fe X] A6374/[0 III] A5007 vs. [Fe VII] A6087/[O III] A5007), we have compared the model 
results with the observations (Osterbrock 1977, 1985; Koski 1978; Osterbrock & Pogge 
1985; Shuder & Osterbrock 1981; Murayama 1998) and found that; 1) the most probable 
values of \ogU and a are —2 and —1, respectively, 2) nn ~ 10''^'^ cm~^, and 3) an unusually 
high iron abundance, up to ten times the solar value while keeping the remaining elements 
the solar values, was required to explain the observed very high [Fe VII] A6087/[O III] A5007 
ratios for the Sis (Murayama 1998). To briefly demonstrate these results, a comparison 
between the observations and the models is shown in Figure 1 which is a diagram of [Fe VII] 
A6087/[O III] A5007 vs. [O I] A6300/[O III] A5007. Note that we adopt logC/ = -2 and 
q; = — 1 for all the models presented here. A serious problem is that a wide span in the iron 
abundance extending up to ten times solar is required to explain the observational data 
points. Although the circumnuclear star-forming regions in Seyfert nuclei show a nitrogen 
overabundance of up to a factor of three at most (e.g., Storchi-Bergmann et al. 1996), it 
is quite unlikely that the NLR in the Seyfert nuclei shows a factor of ten overabundance 
of iron because this overabundance cannot be explained in the terms of ordinary star 
formation history in galaxies. We also mention that single cloud models predict a tendency 
towards smaller [O I]/[0 III] with increasing iron abundance. This may be due to the fact 
that when one increases the iron abundance, the contribution of Fe° and Fe+ to the cooling 
within the partly ionized region increases more than that of O*^. 



2.2. Dual-Component Cloud Models 

As shown in MT98, the observed [Fe VII] /[O III] ratios are systematically higher in 
Sis than in S2s. This suggests that the major contributor to the HINER emission in 
Seyferts has a viewing angle dependence in accordance with current unified models of AGN 
(e.g., Antomicci 1993). It is reasonable to consider that the most probable site of this 
HINER emission corresponds to the inner edges of dusty tori. We therefore proceed to 
construct dual- component models in which the inner surface of a torus is introduced as a 
new ionized-gas component in addition to the traditional NLR component. The observed 
[Fe VII] /[O III] ratios of the Sis are ~ 0.1 on the average while those of the S2s are ~ 0.01 
(MT98). It is hence suggested that the typical NLRs have intensity ratios of [Fe VII] /[O III] 
~ 0.01 and that the observed higher ratios in Sis than in S2s are mostly due to the 
contribution from the HINER torus. 

For the HINER torus component, we calculate photoionization models as follows. The 
single-cloud model suggests that the ionization parameter lies in the range of logU ~ — 1.5 
2. As for the electron density, it is often considered that the inner edges of tori have 
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higher electron densities, e.g., Ue ~ 10^~^ cm~^ (Pier & Voit 1995). This higher density is 
also expected in terms of the locally optimally-emitting cloud (LOG) models (Ferguson et 
al. 1997a). We therefore calculated models covering the following ranges of parameters; 
logU = -2.2, -2.0, -1.8, -1.6, and -1.4, and lognn (cm'^) = 6, 6.5, 7, 7.5, 8, 8.5, and 9. 
In order to increase the [Fe VII] /[O III] ratio by one order of magnitude, the contribution 
to the [Fe VII] emission of the torus component must be very high. Because the largest 
[Fe VII]/[0 III] ratio of the observed data is ~ 0.5, [Fe VIl]/[0 III] of the torus component 
must be greater than 0.5. However, we find that ionization-bounded models cannot explain 
the observed large [Fe VII] /[O III] values by simply increasing electron densities up to 10^ 
cm~^. Further, such very high-density models yield unusually strong [O I] emission with 
respect to [O III]. We therefore assume "truncated" clouds with both large [Fe VII] /[O III] 
ratios and little low-ionization lines for the HINER torus. The calculations were stopped 
at a hydrogen column density when [Fe VII] /[O III] = 1. The results are summarized 
in Table 1. For all the models, we also calculated an [Fe X] A6374/[Fe VII] A6087 ratio. 
Since the average [Fe X]/[Fe VII] ratio for Seyferts is ~ 0.4 (e.g.. Pier & Voit 1995), the 
models with [Fe X]/[Fe VII] » 1 are not appropriate. Also, taking into account the critical 
density of [Fe VII], 3.6x10''' cm~^ (De Robertis & Osterbrock 1984), we adopt the model 
with nu = lO'''-^ cm~^ and logU = -2.0, which gives a ratio [Fe X]/[Fe VII] of 0.8 as a 
representative model for the HINER torus. This model has a characteristic thickness of 
I — Nn/riYL ~ 10^^ '^ cm. Therefore, if we regard that the cloud is matter bounded, its size 
is estimated to be ~ 10^^ cm. On the other hand, if we regard the cloud as an ionization 
bounded one, we have to consider a case that the cloud is so dense and dusty that only its 
surface is photoionized. This is indeed the cloud model considered by Pier & Void (1995). 
Their model gives a typical thickness of the surface layer at UV dust optical depth ra ~ 1, 
^layer ~ 6 X 10^^ Fy"^ T4 Xd^ Cm, whcrc F7 is the ionizing energy flux in units of 10^ erg cm~^ 
s~^, T4 is the gas temperature in units of 10^ K, and Xd is the dust-to-gas ratio relative to 
the Galactic value. This is nearly the same as the characteristic length of our best model. 
Therefore our truncated models can be consistently understood in terms either of matter 
bounded models or of ionization bound models. 

Now we can construct dual-component models combining this torus component model 
with the NLR models with a = -1, logC/ = -2, [Fe] = 1, and lognn (cm~^) = 0-7 
(see section 2.1). When we construct a composite model between the NLR and the torus 
component, we have to take into account the covering factor of the torus. It is known that 
the NLR shows a biconical emission-line morphology with a typical semi-opening angle of 
about 30° (e.g., Pogge 1989; Wilson & Tsvetanov 1994; Schmitt & Kinney 1996). This 
semi-opening angle implies a covering factor of the torus, 1 — An/(47r) ~ 0.9 where Afl is 
the full opening solid angle of the NLR cones in steradian. Therefore, one can infer that 
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the relative contribution of the NLR component may be of the order 10 %. In Figure 2, we 
present the results of the dual-component models. Here the lowest dashed line shows the 
results of the NLR component models with a = —1, logf/ = —2, and [Fe] = 1 as a function 
of n-n from 1 cm~^ to 10^ cm^^. If we allow the contribution from the torus component to 
reach up to ~ 50 % in the Seyferts with very high [Fe VII] /[O III] ratios, we can explain all 
the data points without invoking the unusual iron overabundance. Note that the majority 
of objects can be explained by simply introducing a ~ 10 % contribution from the HINER 
torus. The important quantities characterizing the models are summarized in Table 2. The 
model area shown in Figure 2 fits to the observed data points better than that in Figure 1. 
The dual-component models predict the occurrence of a lower cut-off in the [O I]/[0 III] 
ratio, i.e., no object with log [O l]/[0 III] smaller than —1.8, which is consistent with the 
observations. However, as shown in Figure 1, single cloud models of both low density and 
high iron abundance predict on the contrary the existence of AGNs with log [O l]/[0 III] 
< —1.8. This further supports our contention that the dual- component models studied here 
are more successful than the single cloud models. 

3. CONCLUDING REMARKS 

There has been a debate about the origin of the high-ionization lines. The possible 
origins are summarized as follows: 1) a hot collisionally ionized gas with temperatures of 
Tc ~ 10^ K (Oke & Sargent 1968; Nussbaumer & Osterbrock 1970), 2) a cool gas (Tg ~ 
a few xlO^ K and nn ~ 10^ cm^'^) photoionized by the central non-thermal continuum 
emission (Osterbrock 1969; Nussbaumer & Osterbrock 1970; Grandi 1978), 3) a low-density 
(e.g., ~ 1 cm~^) interstellar medium (ISM) photoionized by the central non-thermal 
continuum emission (Korista & Ferland 1989; see also Murayama, Taniguchi, & Iwasawa 
1998), and 4) a combination of shocks and photoionization by the central non-thermal 
continuum emission (Viegas-Aldrovandi & Gontini 1989). Although photoionization models 
seem to be a modest idea, they cannot explain the observed strong high-ionization lines 
(Murayama, Taniguchi, & Iwasawa 1998 and references therein). However, MT98 found 
that the Sis have excess [Fe VII] A6087 emission with respect to the S2s, implying that the 
majority of the HINER emission in the Sis arises from the inner wall of dusty tori. Indeed 
we have shown that our new dual-component photoionization models incorporating the 
contribution from the HINER torus can explain the observed higher [Fe VII] /[O III] ratios 
in the Sis. Although shock models can also explain the observations (Dopita & Sutherland 
1995, 1996; Dopita et al. 1997; Terlevich et al. 1992), it is important to construct more 
realistic multi-component photoionization models in which detailed physical properties of 
ionized clouds are taken into account (e.g., ionization-bounded and matter-bounded clouds: 
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Binette et al. 1996, 1997; Wilson et al. 1997). 
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Table 1. A summary of the torus models 



lognn 


logU 


log L[oiII] = log I/[FeVII] 


logiVn 


[Fe X]/[Fe VII] ^ 


(cm~^) 




(erg s~^) 


(cm~^) 




7.0 


-1.4 


38.8 


21.1 


53 


7.0 


-1.6 


39.1 


20.8 


12 


7.0 


-1.8 


39.1 


20.3 


3.0 


7.5 


-1.4 


38.9 


21.3 


57 


7.5 


-1.6 


39.2 


21.0 


14 


7.5 


-1.8 


39.4 


20.8 


3.4 


7.5 


-2.0 


39.5 


20.4 


0.8 


R 


—1 4 


38 6 


21 3 


1 1 1 

J. J. J. 


8.0 


-1.6 


39.0 


21.1 


21 


8.0 


-1.8 


39.3 


20.3 


5.2 


8.0 


-2.0 


39.3 


20.6 


1.5 


8.0 


-2.2 


39.1 


20.1 


0.4 


8.5 


-1.4 


38.3 


21.3 


163 


8.5 


-1.6 


38.6 


21.1 


52 


8.5 


-1.8 


38.9 


20.9 


10 


8.5 


-2.0 


39.0 


20.6 


2.8 


8.5 


-2.2 


38.8 


20.3 


0.7 


9.0 


-1.6 


38.2 


21.1 


77 


9.0 


-1.8 


38.4 


20.9 


20 


9.0 


-2.0 


38.5 


20.6 


5.3 


9.0 


-2.2 


38.4 


20.3 


1.4 



^The intensity ratio of [FeX]A6374/[FeVII]A6087. 
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Table 2. Results of the dual-component models 



/[oiii] (Torus) 
/[Oiii] (Total) 


log L[oiii] 
(erg s"^) 


logIv[oi] 

(erg s~^) 


log L[FcVII] 

(erg s~"^) 


lo^r 1 


lo^rt^"^"^ 

[OIIIJ > 


Model 


NLR 




42.3 


41.1 


40.1 


-1.2 


-2.2 


a 


Torus 




39.5 




39.5 







1) 


NLR + Torus 


0.02 


41.2 


40.0 


39.6 


-1.2 


-1.6 


c 


0.10 


40.5 


39.2 


39.5 


-1.3 


-1.0 


d 


0.20 


40.2 


38.9 


39.5 


-1.3 


-0.7 


e 


0.50 


39.8 


38.3 


39.5 


-1.5 


-0.3 


f 



4ognH (cm ^) = 4.0 and logU = —2.0. 

•^logriH (cm-=^) = 7.5, logC/ = -2.0, and logA^n (cm-^) = 20.4. 
'^NLR X 0.078 + Torus. 
^NLR X 0.014 + Torus. 
^NLR X 0.0064 + Torus. 
^NLR X 0.0016 + Torus. 
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Fig. 1. — Single-cloud photoionization models are shown in this diagram of 
[OI]A6300/[OIII]A5007 vs [FcVII]A6087/[OIII]A5007. The circles arc Sis, the triangles are 
SI. 5s, and the squares are S2s. The filled symbols are the objects with [FeX] emission while 
the open symbols are the objects without [FeX] emission. The numbers labeling the lowest 
dashed line represent lognn- 

Fig. 2. — New dual-component photoionization models are shown in this diagram of 
[OI]A6300/[OIII]A5007 vs [FeVII]A6087/[OIII]A5007. Symbols have the same meaning as 
in Figure 1. The numbers labeling the lowest dashed line represent lognn- 



